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ABSTRACT 



Aims. We re-analyse the HARPS radial velocities of HD 10180 and calculate the probabilities of models with differing numbers of 
periodic signals in the data. We test the significance of the seven signals, corresponding to seven exoplanets orbiting the star, in the 
Bayesian framework and perform comparisons of models with up to nine periodicities. 

Methods. We use posterior samplings and Bayesian model probabilities in our analyses together with suitable prior probability 
densities and prior model probabilities to extract all the significant signals from the data and to receive reliable uncertainties for the 
orbital parameters of the six, possibly seven, known exoplanets in the system. 

Results. According to our results, there is evidence for up to nine planets orbiting HD 10180, which would make this this star a record 
holder in having more planets in its orbits than there are in the Solar system. We revise the uncertainties of the previously reported six 
planets in the system, verify the existence of the seventh signal, and announce the detection of two additional statistically significant 
signals in the data. If of planetary origin, these two additional signals would correspond to planets with minimum masses of S.l+H 



and l.9tH M e 



on orbits with 67.55+?™ and 9.655^2™ days periods (denoted using the 99% credibility intervals), respectively. 
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1. Introduction 

Over the recent years, radial velocity surveys or nearby stars 
have provided detections of several exoplanet systems with 
multiple l ow-mass pl a nets, even few Earth-masses, in their or- 
bits (e.g. lLovis et all l2011t iMavor et all I2009L I2011I) . These 
syste ms include a four p lanet system around the M-dwarf GJ 
581 (Ma yor et all 120091) . which has be en proposed to possi- 
bly have a habitable pla net in its orbit (Ivon Paris et all 1201 it 
IWordsworth et all 1201 Ol) . a system of likely a s many as seven 
planets orbiting HD 10180 (lLovis et all 1201 ll) . and several sys- 
tems with 3-4 low mass planets, e.g. HD 20 792 with minimum 
planetary masses of 2.7, 2.4, and 4.8 M e dPepe et all 1201 ll) 
and HD 69830 w ith three Neptune mass planets in its orbit 
(lLovis et al l l2006h . 

Currently, one of the most accurate spectrographs used in 
these surveys, is the High Accuracy Radial Velocity Planet 
Searcher ( HARPS) mounted o n the ESO 3.6m telescope at La 
Silla, Chile dMavor et alll2003"h . In this article, we re -analyse the 
HARP S radial velocities of HD 10180 published in lLovis et aD 
(2011). These measurements were reported to contain 6 strong 
signatures of low-mass exoplanets in orbits ranging from 5 days 
to roughly 2000 days and a possible seventh signal at 1.18 days. 
These planets include five 12 to 25 M ffi planets classified in the 
category of Neptune-like planets, a more massive outer planet 
with a minimum mass of 65 M®, and a likely super-Earth with 
a minimu m mass of 1 .35 M ffi orbiting the star in close proxim- 
ity (lLovis et all 1201 ll) . While the confidence in the existence of 
the six more massive companions in this system is rather h igh, it 
is less so for the innermost super- Earth (Fe roz et all 1201 lb . Yet, 
even if the radial velocity signal corresponding to this low-mass 
companion was an artefact caused by noise and data sampling 
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or periodic phenomena of the stellar surface, the HD 10180 sys- 
tem would be second only to the Solar system with respect to 
the number of planets in it s orbits, together wit h the transiting 
Kepler- 1 1 6-planet system (Lissauer et all 1201 ll) . 

In this article, we re-analyse the radial velocity data of HD 
10180 using posterior samplings and mode l probabilities. We 
perform these analyses to verify the results of lLovis et al. l (l20lTl) 
with another data analysis method, to calculate accurate uncer- 
tainty estimates for the planetary parameters, and to see if this 
data set contains additional statistically significant periodic sig- 
nals that could be interpreted as being of planetary origin. 



2. Observations of HD 10180 planetary system 

The Gl V star HD 10180 is a relatively nearby and bright tar- 
get for radial velocity surveys with a Hipparcos parallax of 
25.39+0.62 mas and V = 7.33 (lLovis et all r201 ll) . It is a very 
inactive (log/? hk - -5.00) Solar-type star with similar mass and 
metallicity (m* = 1.06 ± 0.05, [Fe/H]= 0.08 ± 0.01) and does 
not appear to show an y well-defined act ivity cycles based on the 
HARPS observations (lLovis et alll201 ll) . When ann ouncing the 
discov ery of the planetary system around HD 10180, lL~ovis et al.l 
(1201 ll) estimated the excess variations in the HARPS radial ve- 
locities, usually referred to as stellar jitter, to be very low, ap- 
proximately 1.0 ms _I . These properties make this star a suitable 
target for radial velocity surveys and enable the detection of very 
lo w-mass planets in it s orbit. 

lLovis et all (1201 ll) announced in 2010 that HD 10180 is host 
to six Neptune-mass planets in its orbit with orbital periods of 
5.76, 16.36, 49.7, 123, 601, and 2200 days, respectively. In ad- 
dition, they reported a 1 . 1 8 days power in the periodogram of the 
HARPS radial velocities that, if caused by a planet orbiting the 
star, would correspond to a minimum mass of only 35% more 
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than that of the Earth. These claims were based on 190 HARPS 
measurements of the variations in the stellar radial velocity be- 
tween November 2003 and June 2009. 

The HARPS radial vel ocities have a base line of more than 
2400 days, which enabled lLovis et all d201 lb to constrain the 
orbital parameters of the outer companion in the system with al- 
most similar orbital period. In addition, these HARPS velocities 
have an estimated average instrument uncertainty of 0.57 ms~' 
and a relatively good phase coverage with only seven gaps of 
more than 100 days, corresponding to the annual visibility cycle 
of the star in Chile. 



3. Statistical analyses 

We analyse the HD 10180 radial velocities using a simple model 
that contains k Keplerian signals assumed to be caused by non- 
interacting planets orbiting the star. We also assume that any 
post-Newtonian effects are negligible in the timescale of the ob- 
serva tions. Our statistical mod els are then those described in 
e.g. iTuomi & Kotirantal d2009l) and iTuomil d201 If) , where each 
radial velocity measurement was assumed to be caused by the 
Keplerian signals, some unknown reference velocity about the 
data mean, and two Gaussian random variables with zero means 
representing the instrument noise with a known varian ce as re- 
ported for the HD 10180 data by lLovis et all (1201 lh together 
with the radial velocities, and an additional random variable with 
unknown variance that we treat as a free parameter of our model. 
This additional random variable describes the unknown excess 
noise in the data caused by the instrument and the telescope, at- 
mospheric effects, and the stellar surface phenomena. 

Clearly, the assumption that measurement noise has a 
Gaussian distribution might be limiting in case it was actually 
more centrally concentrated, had longer tails, was skewed, or 
was dependent on time and other possible variables, such as 
stellar activity levels. However, with "only" 190 radial veloci- 
ties it is unlikely that we could spot non-Gaussian features in the 
data reliably. For this reason, and because as far as we know the 
Gaussian one is the only noise model used when analysing radial 
velocity data, we restrict our statistical models to Gaussian ones. 



3.1. Posterior samplings 

We analyse the radial velocities of HD 1018 using the adaptive 
Metropolis posterior sampling algorithm of lHaario et al.l (1200 lh 
because it appears to be efficient in drawing a statistically rep- 
rese ntative sample from the parameter posterior density in prac- 
tice dTuomiluOl lUTuomi et al.L 1201 lb . This algorithm is essen- 
tialy an a daptive version of the famous Metropo lis-Hastings al- 
gorithm (Metropolis et al.Lll953tlHastin^sL IT970). that adapts the 
proposal density to the information gathered up to the ith mem- 
ber of the chain when proposing the i + 1th member. It uses a 
Gaussian multivariate proposal density for the parameter vector 
6 and updates its co variance matrix C,- + i using 



Ci+l = — ; — Ci + - 
I I 



+ el 



(1) 



where 6 is the mean of the parameter vector, (-) r is used to de- 
note the transpose, / is identity matrix of suitable dimension, e is 
some very small number that enables the correct ergodicity prop- 
erties of the resulting chain (Haario et al., 2001), and parameter 
s is a scaling parameter that can be chosen as 2A 2 K~ l , where K 
is the dimension of 9, to optimise the mixing properties of the 
chain (lGelmanetal.Lll996h . 



We calculate the marginal integrals needed in model se- 
lection using the samples from posterior probability densities 
with the one blo c k Me tropolis-Hastings (OB MH) me t hod o f 
IChib & Jeliazkovl (l200lh . also discussed in lClvde et all d2007h . 
However, since the adaptive Metropolis algorith m is not ex- 
actly a Markovian pro cess, only asymptotically so (lHaario et al.L 
l200 lh . the method of IChib & Jeliazkovl (l200lh does not neces- 
sarily yield reliable results. Therefore, after a suitable burn-in 
period used to find the global maximum of the posterior density, 
during which the proposal density also converges to a multivari- 
ate Gaussian that approximates the posterior, we fix the covari- 
ance matrix to its present value, and continue the sampling with 
the Metropolis-Hastings algorithm, that enables the applicability 
of the OBMH method. 



3.2. Prior probability densities 

The prior probability densities of Keplerian models describing 
radial velocity data hav e received little attention in the literature. 
iFord & Gr egory (2007) proposed choosing the Jeffreys' prior for 
the period (P) of the planetary orbit, the radial velocity am- 
plitude (K), and the amplitude of "jitter", i.e. the excess noise 
in the measurements (07). This choice was justified because 
they make the log arithms of thes e param eters evenly distributed 
(Ford & GregorvT l2"007l: iGregorvl E007ab . We use this functional 
form of prior densities for the orbital period and choose the cut- 
off periods such that they correspond to the 1 day period, below 
which we do not expect to find any signals in this work, and a pe- 
riod of 10T o b s , where T t, s is the time span of the observations. 
We choose this upper limit because it enables the detection of 
linear trends in the data corresponding to long-period compan- 
ions whose orbital period cannot be constrai ned, but is no t much 
greater than necessary in practice (see e.g. ITuomi et al.l [2009, 
for the detectability limits as a function of T t, s ), which could 
slow down the posterior samplings by increasing the hypervol- 
ume of the parameter space with reasonably high likelihood val- 
ues. Anyhow, if the period of the outermost companion cannot 
be constrained from above, it would violate our detection crite- 
rion of the p revious subsection. 

Unlike in lFord & Gre"gorvl d2007b . we do not use the Jeffreys' 
prior for the radial velocity amplitudes nor the excess noise pa- 
rameter. Instead, because the HARPS data of HD 10180 deviate 
about their mean less than 20 ms _1 , we use uniform priors for 
these parameters as n(K,) = n{crj) - U(Q,aRv), for all i, and 
use a similar uniform prior for the reference velocity (y) about 
the mean as n(y) = U(-cirv, arv), where we choose the param- 
eter of these priors as a^v = 20 ms~'. While the radial velocity 
amplitudes could in principle have values greater than 20 ms~' 
while the data would still not deviate more than that about the 
mean, we do not consider that p ossibility a feasible one. 

Following IFord & Gregory! (120071) . we choose uniform pri- 
ors for the two angular parameters in the Keplerian model, 
the longitude of pericentre (<w) and the mean anomaly (Mo). 
However, we do not set the prior of orbital eccentricity (e) equal 
to a uniform one between and 1 . Instead, we expect high ec- 
centricities to be less likely in this case because there are at 
least six, likely as many as seven, known planets ofbiting HD 
10180. Therefore, high eccentricities would result in instabil- 
ity and therefore we do not consider their prior probabilities 
to be equal to the low eccentricity orbits. Our choice is then a 
Gaussian prior for the eccentricity, defined as 7r(e,) oc A/ r (0, cr 2 e ) 
(with the corresponding scaling in the unit interval), where the 
parameter of this prior model is set as <r e = 0.3. This choice 
penalises the high eccentricity orbits in practice but still enables 
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them if the data insists so. In practice, with respect to the weight 
this prior puts on zero eccentricity, it gives the eccentricities of 
0.2, 0.4, and 0.6 relative weights of 0.80, 0.41, and 0.14, indicat- 
ing that this prior can only have a relatively minor effect on the 
posterior densities. 

Finally, we required that the planetary systems correspond- 
ing to out Keplerian solution to the data did not have orbital 
crossings between any of the companions. We used this condi- 
tion as additional prior information by estimating that the like- 
lihood of having any two planets in orbits that cross oneanother 
is zero. We could have used a more restrictive criteria, such as 
the requirement that the planets do not enter each others Hill 
spheres at any given time, but decided to keep the situation sim- 
pler because we wanted to see whether the orbital periods of the 
proposed companions get constrained by data instead of stability 
criteria, as described in the next subsection. 

This choice of restricting the solutions in such a way that 
the corresponding planetary system does not suffer from desta- 
bilising orbital crossings also helps reducing the computational 
requirements by making the posterior samplings simpler. Having 
found k Keplerian signals in the data, we simply search for ad- 
ditional signals between them by limiting the period space of 
the additional signals between these k periods. We set the initial 
periods of the k planets close to the solution of the fe-Keplerian 
model and perform k + 1 samplings where each begins with the 
period of the k + 1th signal in different "gaps" around the pre- 
viously found k signals, i.e. corresponding to planet inside the 
orbit of the innermost one, between the two innermost ones, 
and so forth. If a significant k + 1th signal is not found in one 
of these "gaps", the corresponding solution can simply be ne- 
glected. However, if there are signals in two or more gaps, it 
is straightforward to determine the most significant one because 
they can be treated as different models containing the same exact 
number of parameters. We then choose the most significant peri- 
odicity as the k + 1th one and continue testing whether there are 
additional signals in the data. In this way, the problem of being 
able to rearrange the signals in any order, that wou ld cause the 
poste rior density to be actually highly multimodal dFeroz et al.L 
1201 ll) . actually disappears because in a given solution the orbital 
crossings are forbidden and the ordering of the companions re- 
mains fixed. 

3.3. Detection threshold 

While the Bayesian model probabilities can be used reliably in 
assessing the relative posterior probabili ties of models with dif- 
fering numbers of Keplerian s i gnals (e.g. Ford & Gregory, 2007; 
Gregoryl l20lTl lLoredo eTail 1201 ll iTuomi & Kotirantal. 12009: 
TuomiLl201 ll) . we introduce additional criteria to make sure that 
the signals we detect can be interpreted as being of planetary 
origin and not arising from unmodelled features in the measure- 
ment noise or as spurious signals caused by measurement sam- 
pling. Our basic criterion is that the posterior probability of a 
model with k + 1 Keplerian signals has to exceed 150 times that 
of a model with only k Keplerian signals to claim there are k + 1 
planets orbiting the target star. We c hoose this threshold p roba- 
bility based on the considerations of lKass & Raft erv ( 1995). 

We require that the signals we detect in the measurements 
have radial velocity amplitudes, Ki for all i, statistically distin- 
guishable from zero. In practice this means, that not only the 
maximum a posteriori estimate is clearly greater than zero, but 
that t he corresponding Bayesi an 5 credibility sets, as defined in 
e.g. lTuomi & Kotirantal d2009l) . do not allow the amplitude to be 
negligible with 5 = 0.99, i.e. with a probability of 99%. The 



second criterion is that the periods of all signals (Pf) are well 
defined by the posterior samples in the sense, that they can be 
constrained from above and below and do not get constrained 
purely by the condition that orbital crossings corresponding to 
the planetary orbits are not allowed. 

To further increase the confidence of our solutions, we do 
not set the prior probabilities of the different models equal in 
our analyses. We suspect a priori, that detecting k + 1 planets 
would be less likely than detecting k planets in any given system. 
In other words, we estimate that any set of radial velocity data 
would be more likely to contain k Keplerian signals than k + 1 . 
Therefore, we set the a priori probabilities of models with k and 
k + 1 Keplerian signals such that P(k) = 2P(k +1) for all k, i.e. 
we penalise the model with one additional planet by a factor of 
two. Because of this subjective choice, if the model with k + 1 
Keplerian signals receives a posterior probability that exceeds 
our detection threshold of being 150 times greater than that of 
the corresponding model with k Keplerian signals, we are likely 
underestimating the confidence level of the model with k + 1 
Keplerian signals relative to a uniform prior. 

Physically, this choice of prior probabilities for different 
models corresponds to the fact that the more planets there are 
orbiting a star, the less stable orbits there will be left. Therefore, 
we estimate that if k planets are being detected, there is naturally 
"less room" for an additional k + 1th companion. However, this 
might be true statistically, not in an individual case, which leaves 
room for discussion. Yet, this and the benefit that we underesti- 
mate the significance of any detected signals encourages us to 
use this prior. 

3.4. Frequentist and Bayesian detection thresholds 

In addition to the Bayesian analyses described in the previous 
subsections, we analyse the r esiduals of each model usin g the 
L omb-Scargle periodograms Combl H976t IScarglel fl982b . As 
in lLovis et all (120111) . we plot the 10%, 1%, and 0.1% false 
alarm probabilities (FAP's) to the periodograms to se e the signif- 
icance of the powers they contain. However, because lLovis et all 
(201 1) calculated the FAP's in a frequentist manner by perform- 
ing random permutations to the residuals while keeping the ex- 
act epochs of the data fixed and by seeing how often this random 
permutation produces the observed powers, we first put this pe- 
riodogram approach to its philosophical context. 

Generating N random permutations of the residuals for each 
model aims at simulating a situation, where it would be pos- 
sible to receive N independent sets of measurements from the 
system of interest and seeing how often the measurement noise 
generates the signals corresponding to the highest powers in 
the periodogram. While they would not have been independent 
even in the case this process tries to simulate, because the exact 
epochs are fixed making the measurements actually dependent 
through the dimension of time (the measurement is actually a 
vector of two numbers, radial velocity and time), this approach 
suffers from another more significant flaw. The uncertainties of 
the signals removed from the data cannot be taken into account, 
which means that the method assumes the removed signals were 
known correctly. Obviously this is not the case even with the 
strongest signals, and even less so for the weaker ones, making 
the process prone to biases. Therefore, while likely producing 
reliable results when the signals are clear and their periods get 
accurate constraints, this method cannot be expected to provide 
reliable results in the case of extremely wea k signals wi t h larg e 
(and unknown) uncertainties. As noted by lLovis et all d201 lb . 
when testing the significance of the 600 days signal, they could 
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not take into account the uncertainties of the parameters of each 
Keplerian signal, that of the reference velocity, nor the uncer- 
tainty in estimating the excess noise in the data correctly. 

The above "frequentist" way of performing the analyses and 
intepreting the consequent results is different from the Bayesian 
one. Because we only received one set of data, we have to base 
all our results on that and not some hypothetical data that would 
have corresponded to repetition of the original measurements. 
The Bayesian philosophy is to inference all the information from 
the data by combining it with our prior beliefs on what might be 
producing it. For instance, as described above when discussing 
our choice of priors, we can expect tightly packed multiplanet 
systems to be more likely to contain planets on close-circular 
orbits than on very eccentric ones. Also, with the powerful pos- 
terior sampling algorithms available, it is possible to take the 
uncertainties in every parameter into account simultaneously, 
which enables the de te ction of weak signals in the data (e.g. 
iGregorvl 120051 l2007albt iTuomi & Kotirantal 120091) and prevents 
the detection of false positives, as happ ened in the case of Gliese 
581 (IVoet et dll20TalCh-egorvll2^TllTuomill20Tlt . 

Yet, despite the above problems in the traditional peri- 
odogram analyses, we take advantage of the power spectra of 
the residuals in our posterior samplings. The highest peaks 
in the periodogram can be used very efficiently together with 
Bayesian methods by using the corresponding periodicities as 
initial states of the Markov chains in the adaptive Metropolis 
algorithm. Because of this choice, the initial parameter vector 
of the Markov chain starts very close to the likely maximum a 
posteriori (MAP) solution, which makes its convergence to the 
posterior density reasonably rapid and helps reducing computa- 
tional requirements. 

4. Results 

When drawing a sample from the parameter posterior density 
and using it to calculate the corresponding model probabilities, 
it became crucial that this sample was a statistically representa- 
tive one. While posterior samplings generally provide a global 
solutions, it is always possible that the chain converges to a lo- 
cal maximum and stays in its vicinity within a sample of finite 
size. To make sure that we indeed received the global solution, 
we calculated several Markov chains starting from the vicinity 
of the apparent MAP solution and compared them to see that 
they indeed corresponded to the same posterior probability den- 
sity. In practice, sampling the parameter spaces was computa- 
tionally demanding because the probability that the parameter 
vectors drawn from the Gaussian proposal density are close to 
the posterior maximum decreases rapidly when the number of 
parameters with non-Gaussian probability densities increases. 
Therefore, while models with 0-6 planets were reasonably easy 
to sample and we received acceptance rates of 0.1-0.3, these 
rates decreased when increasing the number of signals in the 
model further. As a result, for models with 8-9 Keplerian sig- 
nals, the acceptance rates decreased below 0. 1 and forced us to 
increase the chain lenghts by two orders of magnitude from a 
typical 1.0 x 10 6 to as high as 1.0 x 10 8 . 

In the following subsections, the results are based on sev- 
eral samplings that yielded the same posterior densities, and also 
consistent model probabilities. 

4. 1 . The number of significant periodicities 

The posterior probabilities of the different models provide in- 
formation on the number of Keplerian signals (k) in the data. 



Table 1. The relative posterior probabilities of models with 
k = 0, 9 Keplerian signals (vHO given radial velocities of HD 
10180 (or data, d) together with the periods (P s ) of the signals 
added to the model when increasing the number of signals in 
the model by one. Also shown are the logarithmic Bayesian ev- 
idences (P(d\Mk)) and their uncertainties as standard deviations 
and the root mean square (RMS) values of the residuals for each 
model. 



k 


P(M k \d) 


log P(d\M k ) 


RMS [imT 1 ] 


P s [days] 





l.SxlO" 125 


-621.29 ±0.03 


6.29 




1 


1.6xl0-" 4 


-595.15 ±0.04 


5.39 


5.76 


2 


1.5xl0- 98 


-556.87 ± 0.02 


4.34 


123 


3 


4.1X10" 86 


-528.44 ± 0.07 


3.62 


2200 


4 


2.7xl0~ 53 


-452.17 ±0.13 


2.43 


49.8 


5 


7.2xl0~ 20 


-374.42 ± 0.05 


1.59 


16.35 


6 


3.4xl0~ 13 


-358.36 ±0.12 


1.41 


600 


7 


3.8xl0~ 7 


-343.73 ± 0.09 


1.32 


1.18 


8 


0.003 


-334.06 ± 0.06 


1.24 


67.5 


9 


0.997 


-327.79 ± 0.06 


1.18 


9.65 



lLovis et al.l (1201 lh were confident that there are six planets or- 
biting HD 10180 based on their periodogram-based analyses of 
model residuals and the corresponding random permutations of 
them when calculating the significance levels of the periodogram 
powers. They also concluded that the six planets in the system 
with increasing periods of 5.75962±0.00028, 16.3567+0.0043, 
49.747 ±0.024, 122.72±0.20, 602+11, and 2248+J^ days com- 
prise a stable system given that the masses of the planets are 
within a factor of ~ 3 from the minimum masses of 13.70±0.63, 
11.94±0.75, 25.4+1.4, 23.6±1.7, 21.4±3.0, and 65.3±4.6 M e , 
respectively. 

Because ILovis et af] (l201ll) pointed out that there are actu- 
ally two peaks in the periodogram corresponding to the seventh 
signal, i.e. 1.18 and 6.51 days periodicities that are the one day 
aliases of each other, but noted that the 6.51 signal, if corre- 
sponding to a planet, would cause the system to be unstable in 
short timescales, we adopted the 1.18 periodicity as the seventh 
signal in the data. The relative probabilities of the models with 
k = 0, 9 are shown in Table [T]together with the period (P s ) of 
the next Keplerian signal added to the model. 

According to the model probabilities in Table [T] the 8- and 
9-Keplerian models are the most probable descriptions of the 
processes producing the data out of those considered. Improving 
the statistical model by adding the seventh signal, with a pe- 
riod of 1.18 days, increases the model probability by a factor of 
more than 10 6 , which makes the credibility of this seventh signal 
high. Adding two more signals corresponding to 67.6 and 9.66 
days periodicities increases the model probabilities even further. 
As a result, the 9-Keplerian model receives the greatest poste- 
rior probability of slightly more than 150 times more than the 
next best model, the 8-Keplerian model. This enables us to con- 
clude, that there is strong evidence in favour of the 67.6 and 9.66 
days periodicities not being produced by random processes, i.e. 
measurement noise. 



4.2. Periodograms of residuals 

We subtracted the models with 6 to 9 periodic signal s from 
the data and calculat ed the Lomb-Scargle periodograms (Lomb, 
1 19761: IScargleLI 1982b of these residuals (Pig. [TJ t ogether with th e 
standard analytica FAP's. As already seen in lLovis et al. l d20Tll) . 
the two strong powers corresponding to 1.18 days periodicity 
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1 10 100 1000 10000 

Period [d] 

Fig. 1. The Lomb-Scargle periodograms of the HD 10180 radial 
velocities for the residuals of the models with 6 (top) to 9 (bot- 
tom) periodic signals. The dotted, dashed, and dot-dashed lines 
indicate the 10%, 1%, and 0.1% FAPs, respectively. 



and its 1 day alias at 6.51 days are strong in the residuals of 
the 6-Keplerian model (top panel in Fig.Q]) and exceed the 1% 
FAR These peaks get also removed from the residuals of the 7- 
Keplerian model (2nd panel). However, it can be seen that 9.66 
and 67.6 days periods have strong powers in these residuals, yet, 
neither of them exceeds even the 10% FAP level. Modelling 
these periodicities as Keplerian signals and plotting the pe- 
riodograms of the corresponding residuals of the 9-Keplerian 
model (bottom panel) shows that there are no strong powers left 
in the residuals. While this does not indicate that these two peri- 
odicities are significant, it shows that they are clearly present in 
the periodograms of the residuals and support the findings in the 
previous subsection based on the model probabilities. 

We note that there appear to be two almost equally strong 
peaks in the 8-Keplerian model residuals (Fig. Q] panel 3). 
However, these powers corresponding to periods of 9.66 and 
1.11 days are one day aliases of oneanother. This aliasing is clear 
as the 1.11 days power is absent in the periodogram of the 9- 
Keplerian model (Fig.Q~| bottom panel). 



4.3. Planetary interpretation: orbital parameters 

Because of the samplings of parameter posterior densities of 
each statistical model, we were able to calculate the estimated 
shapes of the parameter distributions for each model and use 
these to estimate the features in the corresponding densities. We 
describe these densities using three numbers, the MAP estimates 
and the correspo nding 99% Bayesian credi bility sets (BCS's), 
as defined in e.g. iTuomi & Kotiranta (120091) . The simple MAP 
point estimates and the corresponding 99% BCS's do not rep- 
resent these dentities very well because some of the parameters 
are highly skewed and have tails on one or both sides. However, 
we list the parameters interpreted as being of planetary origin in 
this subsection. 

When calculating the semimajor axes and minimum masses 
of the planets, we took the uncertainty in the stellar mass into 
account by treating it as an independent random variable. We 
assumed that this random variable h ad a Gaussian de nsity with 
mean equal to the estimate given bv lLovis et all J20lTh of 1.06 
M Q and a standard deviation of 5% of this estimate. As a conse- 
quence, the densities of these parameters are broader than they 
would be if using a fixed value for the stellar mass, which indi- 
cates greater uncertainty in their values. 

4.3.1 . The 6 planet solution 

The parameter estimates of our 6-Keplerian model are listed in 
Table [H This solut ion is consistent with the solution reported 
by lLovis et all d201 ll) but the uncertainties are slightly greater 
likely because we took in to account the unc ertainty in the jitter 
parameter crj and because lLovis et all (1201 ll) used more conser- 
vative uncertainty estimates from the covariance matrix of the 
parameters that does not take the nonlinear correlations between 
the parameters into account. The greatest difference is therefore 
in the uncertainty of the orbital period of the outermost com- 
panion, whose probability density has a long tail towards longer 
periods and periods as high as 2670 days cannot be ruled out 
with 99% confidence (the supremum of the 99% BCS). 

4.3.2. The 9 planet solution 

Assuming that all the periodic signals in the data are indeed 
caused by planetary companions orbiting the star, the param- 
eters of out 9-Keplerian solution are listed in Table [3] and the 
phase folded orbits of the 9 Keplerian signals are plotted in Fig. 

m 

Because the simple estimates in Table [3] can be very mis- 
leading in practice, especially in the case of nonlinear correla- 
tions between the parameters, we also plotted the projected dis- 
tributions of some of the parameters in the Appendix. The dis- 
tributions of Pj, Kj, and e,, for each planet i = 1, ...,9 (with in- 
dice 1 (9) referring to the shortest (longest) period), are shown 
in Fig. I A. 1 1 and show that the periods of all Keplerian signals 
are indeed well constrained, the radial velocity amplitudes differ 
significantly from zero, and the orbital eccentricities peak at or 
close to zero indicating likely circular orbits. We also plotted the 
parameter describing the magnitude of the excess noise, or jitter, 
in Fig lA.ll to demonstrate that the MAP estimate of this param- 
eter of 1.15 ms _I , with a BCS of [0.92, 1.42] ms~', is consistent 



1 lLovis et al] d2011h use letter b for the 1.18 days periodicity not 
present in the 6-Keplerian model. Therefore, we use letters c - h in 
T able 2 to have the letters denote the same signals. Yet, the solution 
of lFeroz et al] d201ll) denotes the 5.76 days signal with letter b. 
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Table 2. The 6-planet solution of HD 10180 radial velocities. MAP estimates of the parameters and their 99% BCS's. 



Parameter 


HD 10180 c 


HD 10180 d 


HD 10180 e 


P [days] 


5.7596 [5.7588, 5.7606] 


16.354 [16.340, 16.368] 


49.74 [49.67, 49.82] 


e 


0.06 [0.00, 0.17] 


0.12 [0.00, 0.29] 


0.01 [0.00, 0.14] 


K [ms- 1 ] 


4.54 [4.06, 5.02] 


2.85 [2.35, 3.34] 


4.32 [3.74, 4.83] 


a) [rad] 


5.8 [-] 


5.8 [-] 


2.6 [-] 


Mo [rad] 


4.7 [-] 


5.8 [-] 


3.7 [-] 


m p sin i [M e ] 


13.2 [11.2, 15.1] 


11.8 [9.5, 14.2] 


25.6 [21.5, 29.7] 


a [AU] 


0.0641 [0.0608, 0.0673] 


0.1287 [0.1220, 0.1349] 


0.270 [0.256, 0.283] 




HD 10180 f 


HD 10180 g 


HD 10180 h 


P [days] 


122.76 [122.152, 123.44] 


603 [568, 642] 


2270 [2020, 2670] 




n 06 rn nn n 

U.UU [U.UU, U.Z.UJ 


n rn nn f) zLQl 


n rn rn nn n t, 1 i 


K [ms- 1 ] 


2.88 [2.28, 3.42] 


1.46 [0.71,2.21] 


2.97 [2.40, 3.68] 


to [rad] 


5.8 [-] 


6.0 [-] 


2.8 [-] 


M [rad] 


4.9 [-] 


4.5 [-] 


3.7 [-] 


m p sin ( [M e ] 


23.1 [18.2, 28.4] 


19.4 [10.2, 29.8] 


64.5 [51.5,78.9] 


a [AU] 


0.491 [0.468, 0.518] 


1.42 [1.33, 1.52] 


3.44 [3.15, 3.88] 


y [ms -1 ] 


-0.43 [-0.87, -0.02] 






a"; [ms -1 ] 


1.40 [1.18, 1.70] 







Table 3. The 9-planet solution of HD 10180 radial velocities. MAP estimates of the parameters and their 99% BCS's. 



Parameter 


HD 10180 b 


HD 10180 c 


HD 10180 i 


P [days] 


1.17766 [1.17744, 1.17787] 


5.75973 [5.75890, 5.76047] 


9.655 [9.583,9.677] 


e 


0.05 [0.00, 0.54] 


0.07 [0.00, 0.16] 


0.05 [0.00, 0.28] 


K [ms -1 ] 


0.78 [0.34, 1.21] 


4.50 [4.07, 4.92] 


0.53 [0.02, 0.99] 


to [rad] 


0.7 [-] 


5.8 [-] 


2.4 [] 


M [rad] 


1.2 [-] 


4.7 [-] 


1.6 [-] 


m p sin i [M e ] 


1.3 [0.5, 2.1] 


13.0 [11.2, 15.0] 


1.9 [0.1,3.5] 


a [AU] 


0.0222 [0.0211, 0.0233] 


0.0641 [0.0610, 0.0671] 


0.0904 [0.0857, 0.0947] 




HD 10180 d 


HD 10180 e 


HD 10180 j 


P [days] 


16.354 [16.341, 16.366] 


49.75 [49.68, 49.82] 


67.55 [66.67, 68.23] 


e 


0.11 [0.00, 0.24] 


0.01 [0.00, 0.11] 


0.07 [0.00, 0.19] 


K [ms -1 ] 


2.90 [2.45, 3.34] 


4.14 [3.67, 4.66] 


0.75 [0.29, 1.21] 


t±> [rad] 


5.6 [-] 


1.6 [-] 


1.6 [-] 


M [rad] 


6.0 [-] 


2.0 [-] 


6.0 [-] 


m p sini [M e ] 


11.9 [9.9, 14.2] 


25.0 [21.1,28.9] 


5.1 [1.9, 8.2] 


a [AU] 


0.1284 [0.1223, 0.1346] 


0.270 [0.257, 0.283] 


0.330 [0.314, 0.347] 




HD 10180 f 


HD 10180 g 


HD 10180 h 


P [days] 


122.88 [122.28, 123.53] 


596 [559, 626] 


2300 [2010, 2850] 


e 


0.13 [0.00, 0.28] 


0.03 [0.00, 0.43] 


0.18 [0.00, 0.44] 


K [ms" 1 ] 


2.86 [2.33, 3.38] 


1.61 [1.00, 2.23] 


3.02 [2.47, 3.64] 


a> [rad] 


5.8 [-] 


0.3 [-] 


4.1 [-] 


M [rad] 


4.9 [-] 


2.6 [-] 


4.1 [-] 


m p sin i [M ffi ] 


22.8 [18.6, 27.9] 


22.0 [13.3, 30.8] 


65.8 [52.9, 78.7] 


a [AU] 


0.494 [0.468, 0.517] 


1.415 [1.324, 1.498] 


3.49 [3.13, 4.09] 


y [ms -1 ] 


-0.34 [-0.75, -0.01] 






o-j [ms -1 ] 


1.15 [0.92, 1.42] 







with the estimate of lLovis et al. of 1 .0 ms , while this is 

not the case for the 6-Keplerian model for which the cr, receives 
a MAP estimate of 1.40 ms" 1 with a BCS of [1.18, 1.70] ms" 1 
(Tabled. 



The periods of all companions get constrained well but that 
of the 9.66 days signal remains bimodal with two peaks at 9.66 
and 9.59 days (Fig. IA.lt . The 9.59 days peak was found to have 
lower probability based on our posterior samplings because the 
Markov chain had good mixing properties in the sense that it 
visited both maxima several times during all the samplings, the 
posterior density had always a global maximum at 9.66 days. 



4.4. Dynamically allowed orbits 



We perform tests of dynamical stability within the context of 
Lagrange stability to see whether the two additional signals in 
the HD 10180 radial velocities could correspond to low-mass 
planets orbiting the star. We u se the analytical app r oxim ated 
Lagrange stability criterion of iBarnes & Greenberg (2006JI to 
test the stability of each subsequent pairs of planets in the sys- 
tem. While this analytical criterion is only a rough approxima- 
tion and only applicable for two planet systems, it can neverthe- 
less provide useful information o n the likely stability or insta - 
bility of the system. According to lBarnes & Greenbergl (12006). 
the orbits of two planets (denoted using subindices 1 and 2, re- 
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Fig. 2. The phase-folded Keplerian signals of the 9-planet solution with the other 8 signals removed. 



spectively) satisfy approximately the Lagrange stability criterion 
if 



where yu,- = m,M 1 , a - + ^ 7i = -^1 - e?, 5 = yjai/a\, 
M - m* + nt\ + m.2, e-, is the eccentricity, a, is the semimajor 
axis, nii is the planetary mass, and ra* is stellar mass. 

Using the above relation, we calculate the threshold curves 
for the z'th planet with both, the next planet inside its orbit (i— 1th) 
and the next planet outside its orbit (/ + 1th). We use the MAP 
parameter estimates for the /th planet and calculate the allowed 
eccentricities of the /—1th and /+ 1th planets as a function of their 
semim ajor axes by using the MAP estimates for their masses. 

As lLovis et all d201 ll) found their 6 planet solution stable, 
we use it as a test case when calculating the Lagrange stability 
threshold curves. We use the 6 planet solution in Table [2] and 
plot the threshold curves together with the orbital parameters in 
Fig. [3] In this Fig., the shaded areas indicate the likely unstable 
parameter space and the red circles indicate the positions of the 
modelled planets in the system. 

It can be seen in Fig. [3] that the /th planet has orbital pa- 
rameters that keeps it inside the Lagrange stability region of the 
neighbouring planetary companions for all / = 1, 6. This re- 
sult then agrees with the numerical integrations of lLovis et alJ 




Fig. 3. The approximated Lagrange stability thresholds between 
each two planets and the MAP orbital parameters of the 6 planet 
solution (Table 13 . 



(1201 ll) and, while only a rough approximation of the reality, en- 
courages us to use the criterion in Eq. (|2} for our 7-, 8-, and 
9-companion solutions as well. Fig. [3] also suggests, that there 
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Fig. 4. The approximated Lagrange stability thresholds between 
each two planets and the MAP orbital parameters of the 9 planet 
solution (Table [3). 



might be stable regions between the orbits of these 6 planets for 
additional low-mass companions. 

When interpreting all the signals in our 9-planet solution as 
being of planetary origin, the stability thresholds show some in- 
teresting features (Fig.|4). The periodicities at 9.66 and 67.6 days 
would correspond to planets that satisfy the condition in Eq. (0 
if the orbital eccentricities of all the companions were close to 
or below the MAP estimate, which, according to the probabil- 
ity densities in Fig. IA.1I appears to be likely based on the data 
alone. This means that the planetary origin of these periodicities 
cannot be ruled out by this analysis. 

In reality, the stability constraints are necessarily more lim- 
iting than those described by the simple Eq. [2] because of the 
gravitational interactions between all the planets, not only the 
nearby ones. Also, it does not take the stabilising or destabil- 
ising effects of mean motion resonances. However, the numeri- 
cal integrat i ons of the orbits of the seven planets performed by 
lLovis et all (1201 lh . actually do not rule the 0.09 and 0.33 AU or- 
bits (Table [3) out as unstable but show that there are regions of 
at least "reasonable stability" in the vicinity of these orbits given 
that they are close-circular and that the planetary masses in these 
orbits are small. When interpreted as being of planetary origin, 
the periodic signals at roughly 9.66 and 67.6 days satisfy these 
requirements. 

We note, that the selected prior density for the orbital eccen- 
tricities, namely 7r(e,) oc 7V(0, cr^) for all i, in fact helps slightly 
in removing a priori unstable solutions from the parameter pos- 
terior density. However, this effect is not very significant in this 
case, because the prior requirement that the signals in the data 
do not correspond to planets with crossing orbits constrains the 
eccentricities much more strongly and the parameter posteriors 
we would receive with uniform eccentricity prior would there- 
fore not differ significantly from those reported in Table [3] and 
Figs.l2landlA~T1 

4.5. Avoiding unconstrained solutions 

To further emphasise our confidence in the 9 periodic signals in 
the data, we tried finding additional signals in the gaps between 
the 9-Keplerians, especially, between the 123 and 600 days or- 
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Fig. 5. The posterior densities of the minimum mass and orbital 
period of the planet that could exist in the habitable zone of HD 
10180 without having been detected using the current data. The 
solid curve is a Gaussian density with the same mean (jx) and 
variance (<x 2 ) as the parameter distribution has. 



bits. This part of the period-space is interesting because any hab- 
itable planet in the system would have its orbital period in this 
space and because the stability thresholds allow the existence of 
low-mass planets in close-circular orbits in this region (see Fig. 
0. 

As could already be suspected based on the periodograms of 
residuals in Fig.[Tj we could not find any signals between the pe- 
riods of 123 and 600 days. The sampling of the parameter space 
of this 10-Keplerian model was much more difficult than that of 
the models with fewer signals because the orbital period of this 
hypothetical 10th signal was only constrained by the fact that a 
priori we did not allow orbital crossings (Fig. |5). As a result, 
the probability density of the orbital period did not have a clear 
maximum but several small maxima, whose relative significance 
is not known because we cannot be sure whether the Markov 
chain converged to the posterior in this case (Fig. |5). Therefore, 
the distribution of the orbital period in Fig. [5] is only a rough 
estimate of what the density might look like. 

Because different samplings yielded similar but not equal 
densities for the orbital period, we could not be sure whether the 
chain had indeed converged to the posterior or not. For this rea- 
son, we did not consider the corresponding posterior probability 
of this model trustworthy and do not show it in Table [TJ The 
posterior probabilities we received were roughly 1-5% of that of 
the 9-Keplerian model. However, these samplings still provide 
some interesting information in the sense that we can put an up- 
per limit to the planetary masses that could exist between the 123 
and 600 days periods and still not be detected confidently by the 
current observations. 
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According to the samplings of the parameter space, the prob- 
ability of there being Keplerian signals between 123 and 600 
days periods with radial velocity amplitudes in excess of 1.1 
ms~' is less than 1%. Our MAP solution for this amplitude is 
0.12 ms" 1 with a BCS of [0.00, 1.10] ms -1 . This means that we 
can rule out the existence of planets more massive than approx- 
imately 12.1 M e in this period space because such companions 
could have been detected by the current data. Therefore, even 
the fact that a signal was not detected can help constraining the 
properties of the system, as seen in the probability density of the 
minimum mass of this hypothetical planet in Fig. [5]- clearly this 
signal is indistinguishable from one with negligible amplitude, 
as the density is peaking close to zero. This result means that an 
Earth mass planet could exist in the habitable zone of HD 10180 
and most likely, if there is a low-mass companion orbiting the 
star between orbital periods of 123 and 600 days, it has a mini- 
mum mass of less than 12.1 M ffi . However, we cannot say much 
about the possible orbit of this hypothetical companion, because 
all the orbits between 123 and 600 days without orbital cross- 
ings are almost equally probable (Fig. [5). All we can say, is that 
orbital crossings limit the allowed periodicities and yield a 99% 
BCS of [128, 534] days for this orbital period, though, further 
dynamical constraints would na rrow this interval e ven further, 
as shown in Fig.|4]and Fig. 12 of lLovis et al. 



4.6. Detectability of Keplerian signals 

To further emphasise the significance of the 9 signals we de- 
tect in the HD 10180 radial velocities, we generated an artificial 
data set to see if known signals could be extracted from it confi- 
dently given the definition of our criteria for detection threshold. 
This data set was genera ted by using the s ame 190 epochs as 
in the HD 10180 data of lLovis et aTTd201 lh . We generated the 
radial velocities corresponding to these epochs by using a super- 
position of the 9 signals with parameters roughly as in Table [3] 
Further, we added three noise components, Gaussian noise with 
zero mean as describe d by the uncertainty estimates of each orig- 
inal radial velocity of lLovis et alj d201 lh . Gaussian noise with 
zero mean and <x = 1.1 ms~' and uniform noise as a random 
number drawn from the interval [-0.2, 0.2] ms -1 . The lattter two 
produce together the observed excess noise in the data of roughly 
1.15 ms -1 when modelled as pure Gaussian noise. We used this 
different noise model when generating the data to not commit an 
inverse crime, i.e. to not generate the data using the same model 
used to analyse it which would co rrespond to studying the prop- 
erties of the model only (see e.g. Kai pio & Somersalol l2005h . 

Analysing this artificial data yielded results confirming the 
trustworthiness of our methods. Using the detection criteria de- 
fined above, we could extract all 9 signals from the artificial data 
with well constrained amplitudes and periods that were consis- 
tent with the added signals in the sense that the 99% BCS's of the 
parameters contained the values of the added signals. The MAP 
estimates did differ from the added signals but not significantly 
so given the uncertainties as described by distributions corre- 
sponding to the parameter posterior density. This simply rep- 
resents the statistical nature of the solutions based on posterior 
samplings. As an example, we show the periods and amplitudes 
of the three weakest signals with the lowes radial velocity am- 
plitudes as probability distributions (Fig. |6). This Fig. indicates 
that if the radial velocity noise is indeed dominated by Gaussian 
noise, the low-amplitude signals we report can be detected confi- 
dently. This conclusion was also supported by the corresponding 
model probabilities we received for the artificial data set. These 
probabilities indicated that the 9-Keplerian model had the great- 
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Fig. 6. The distributions of the periods and amplitudes received 
for the three weakest signals in the artificially generated data. 
The added signals had P, = 1.18, 9.66, and 67.55 days, and K, = 
0.78, 0.53, and 0.75 ms~' for i = 1,3,6, respectively. 



est posterior probability exceedind the threshold of 150 times 
greater than the probability of the 8-Keplerian model. 



4. 7. Comparison with earlier results 



The analyses of lLovis et alJ (1201 ll) of the same radial velocities 
yielded differing results, i.e. the number of significant powers in 
periodogram was found to be 7 instead of the 9 significant peri- 
odic signals reported here. While this difference is likely due 
to the fact that the power spectmms are calculated by fixing 
the parameters of the previous signals to some poin t estimates 
when searching for additional peaks, the analyses of iFeroz et al.1 
d2011l) s uffer from s i milar sources of bias. The Bayesian ap- 
proach of lFeroz et al .1 (1201 lh was basically a search of k+ 1th sig- 
nal in the residuals of the model with k Keplerian signals. They 
derived the probability densities of the residuals by assuming 
they were uncorrelated, which is unlikely to be the case, espe- 
cially, if there are signals left in the residuals. Therefore, any sig- 
nificance test, in this case the comparison of Bayesian evidences 
of a null hypothesis and an alternative one with a model con- 
taining one more signal, is similarly biased by the fact that these 
correlations are not fully accounted for. Whil e this source of bia s 
might be relatively small, in the approach of lFeroz et al 1 (l2011h . 
the effective number of parameters is artificially decreased by 
the very fact that residuals are being analysed. This decrease, 
in turn, might make any comparisons of Bayesian evidences bi- 
ased. Modifying the uncertainties corresponding to the posterior 
density of the model residuals does not necessarily account for 
this decrease in the dimension when the weakest signals among 
the k detected ones are at or below the residual uncertainty and 
their contribution to the total uncertainty of the residuals be- 
comes negligible in the first place. 
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In fact, we could replicate the results of iFeroz et al 1 (120111) 
using the OBMH estimates of marginal integrals. Using their 
simple method, we received a result that the 6-Keplerian model 
residuals could not be found to contain one more signal. While 
the Bayesian evidence for this additional signal did not exceed 
that of having the 6-Keplerian model residuals consist of purely 
Gaussian noise, the amplitude of this additional signal in the 
residuals was also found indistinguishable from zero, in accor- 
dance with our criteria of not detecting a signal. Therefore, the 
results we presen t here do not actually conflict with those of 
Feroz et al. l d20Tlh . As an example, the log-Bayesian evidences 
of 6 and 7-Keplerian models (Table[TJ were -358.36 and -343.73, 
respectively. Analysing the residuals of the model with k = 6 
using 0- and 1 -Keplerian models should yield the si milar num- 
bers ( up to differences in the priors), if the method of lFeroz et al.l 
d201 ll) were trustworthy. Instead, we received values of roughly 
-338 for both models. This means, that the null-hypothesis is 
exaggerated because of the fact that the corresponding model 
contains only two free parameters, the jitter magnitude and the 
reference velocity. The Occamian principle cannot therefore pe- 
nalise this model as much as it should, and the results are biased 
in favour of the null-hypothesis, which effectively prevents the 
detection of low-amplitude signals. 

We also tested the method of Fero z et al.l d201 lb in analysing 
the artificial test data described in the previous subsection. The 
results were almost similar: the log-Bayesian evidences given 
the residuals of a model with k Keplerian signals were found 
to favour the 6-Keplerian interpretation, while the artificial data 
was known to contai n 9 signal s . This further emphasises the fact, 
that the method of lFeroz et al.l (120 1 ll) . while capable of detecting 
the strongest signals in the data, cannot be considered trustwor- 
thy if it fails to make a positive detection of a low-amplitude 
signal. Yet, it is likely trustworthy if it provides a positive detec- 
tion. 

5. Conclusions and Discussion 

We have re-analys ed the 190 HARP S radial velocities of HD 
10180 published in lLovis et al.l d201 ll) and report our findings in 
this article. First, we have revised the orbital parameters of the 
proposed 6 planetary companions to this star and calculated real- 
istic uncertainty estimates based on samplings of the parameter 
space. We a lso verified t he si gnificance of the 1.18 days signal 
reported by lLovis et all (1201 ll) and interpreted as arising from 
a planetary companion with this orbital period and a minimum 
mass of as low as 1.3 M e . In addition to these seven signals, 
we report two additional periodic signals that are, according to 
our model probabilities in Table Q] statistically significant and 
unlikely to be caused by noise or data sampling or poor phase- 
coverage of the observations. Their amplitudes are well con- 
strained and differ statistically from zero, which would not be 
the case unless they corresponded to actual periodicities in the 
data. We can also constrain their periods from above and below 
reasonably accurately. 

A related analysis of t he same radial velocities was recently 
carried out bv lFeroz et al.l (1201 ll) but they received results differ- 
ing in the number of significant periodicities. They claimed that 
only 6 signals can be detected reliably in the data, as opposed 
to 9 detected in the current work. However, they first analysed 
the data using a model with k Keplerian signals and analysed the 
remaining residuals to see if they c ontained one more , k + 1th 
signal. While, as in the analyses of lLovis et al.1 (1201 ll) . this ap- 
proach does not fully account for the uncertainties of the first 
k signals when they have low amplitudes and their contribution 



to the residual uncertainty is negligible, it also actually assumes 
that the fc-Keplerian model is a correct one and then tests if it is 
not so, which is a clear contradiction and, while useful in case 
of strong signals, as demonstrated bv lFeroz et al.l (1201 ll) . likely 
prevents the detection o f weak signals in the data. This is un- 
derlined by the fact that Fer oz et al.l (1201 ll) assume the residual 
vector to have an uncorrected multivariate Gaussian distribu- 
tion - this clearly cannot be the case if there are signals left in 
the residuals. Our analyses are not prone to similar weaknesses. 

Because planetary companions orbiting the star would pro- 
duce the kind of periodicities we observe in the radial velocities, 
the interpretation of the two new signal s as two new low- mass 
planets seems reasonable. As noted by lLovis et all (1201 ll) . the 
star is a very quiet one without clear activity-induced periodici- 
ties, which makes it unlikely that one or some of the periodic sig- 
nals in the data were caused by stellar phenomena. Also, the pe- 
riodicities we report, namely 9.66 and 67.6 days, do not coincide 
with any periodicities arising from the movement of the bodies 
in the Solar system. Therefore, we consider the interpretation 
of these two new signals of being of planetary origin to be the 
most credible explanation. If this was the case, these two signals 
would correspond to planets on close-circular orbits with mini- 
mum masses of 1.9+J'g and M e , respectively, enabling the 
classification of them as super-Earths. 

Apart from the significance of the signals we observe, there 
is another rather strong argument in favour of the interpretation 
that all nine signals in the data are actually of planetary ori- 
gin. Assuming that they were not, which based on stability rea- 
sons is the case with the 6.51 days signal that is quite certainly 
an alias of the 1.1 8 days periodicity likely caused by a planet 
(lLovis et all 1201 ll) . we would expect the weakest signals to be 
at random periods independent of the six strong periodicities in 
the data and the seventh 1.18 periodicity. Instead, this is not the 
case but the two additional signals reported in the current work 
appear at periods that fall in between the existing ones and, if 
interpreted as being of planetary origin, likely have orbits that 
enable long-term stability of the system (Fig. |4]l if their orbital 
eccentrici ties are clo s e to o r below the estimates in Table [3] As 
stated bv lLovis et al.l (1201 ll) . there are "empty" places in the HD 
10180 system that allow dynamical stability of low-mass planets 
in the orbits corresponding to these empty places in the orbital 
parameter space, especially in the a - e space. The two periodic 
signals we observe in the data correspond exactly to those empty 
places if interpreted as being of planetary origin. 

Additional measurements are needed to verify the signifi- 
cance of the two new periodic signals in the radial velocities 
of HD 10180 and to set tighter constraints to the orbital param- 
eters of the planets in the system. Also, the possibility that all 
nine signals in the data correspond to planets should be tested 
by full-scale numerical integrations of their orbits. If all the con- 
figurations allowed by our solution in Table [3] were found to 
correspond to unstable systems in any ti mescales of less th an 
the estimated stellar age of 4.3+0.5 Gyr (lLovis et all 1201 lb . it 
would be a strong argument against the planetary interpretation 
of one or both of the signals we report in this article or the third 
low-amplitude 1. 18 days sign a l. Ho wever, the results we present 
here and those in lLovis et al.l d201 ll) suggest that this planetary 
interpretation of all the signals cannot be ruled out by dynamical 
analyses of the system. 

If the significance of these signals increases when additional 
high-precision radial velocities become available, and their in- 
terpretation as being of planetary origin is confirmed, the plan- 
etary system aroung HD 10180 will be the first one to top 
the Solar system in terms of number of planets in its orbits. 
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Further, according to the rough dynamical considerations of the 
current work and the more extensive numerical integrations of 
lLovis et all d201 ll) . there are stable orbits for a low-mass com- 
panion in or around the habitable zone of the star. If such a com- 
panion exists, its minimum mass is unlikely to exceed 12.1 M ffi 
according to our posterior samplings of the corresponding pa- 
rameter space. 
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Fig. A.l. The posterior densities of orbital periods (Pi), eccentricities (<?;), and radial velocity amplitudes (Ki) corresponding to the 
9 Keplerian signals in the data and of the magnitude of the excess jitter (cry). The solid curve is a Gaussian density with the same 
mean (p) and variance (cr 2 ) as the parameter distribution has. 
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